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r^ I Abstract 

The singular sector of zero genus case for the Hermitian random matrix model in the large N 

limit is analyzed. It is proved that the singular sector of the hodograph solutions for the underlying 

(^ ' dispersionless Toda hierarchy and the singular sector of the 1-layer Benney (classical long wave 

r^ ■ equation) hierarchy are deeply connected. This property is due to the fact that the hodograph 

equations for both hierarchies describe the critical points of solutions of Euler-Poisson-Darboux 

equations E{a,a), with a = —1/2 for the dToda hierarchy and a = 1/2 for the 1-layer Benney 



> 






(~^ ■ hierarchy. 

o' 
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1 Introduction 

The Hermitian matrix (H) model is nowadays the paradigmatic model in the theory of random 
matrices ( see e.g. [1, 2]). The partition function for this model is 

Zn= fdHeKp(^-^tr{Y,ikH')), (1) 

•^ fe>i 

where T > is a real parameter which represents the temperature and the integration is performed 
on the space oi NxN Hermitian matrices. The large N limit of the H model is of particular interest 
since it exhibits many important universality properties. It turns out that most of the applications 
of the H model [1, 2] arise after regularizing the large A^-limit solutions at their singular points 
(double-scaling limit method). In this sense the analysis and characterization of the singular sector 
of the H model is of great relevance. 

The simplest situation corresponds to the so called zero genus case in which the support of 
the eigenvalue density reduces to a single interval [/3i,/32] as A^ — )■ oo. Moreover, it is well-known 
( see e.g. [3, 4]) that the endpoints /3 := (/3i,/32) evolve with the temperature T and the coupling 
constants t^ according to the dispersionless Toda (dToda) hierarchy . The first members of this 
hierarchy are the dToda equation vtt = (logw)iifi or equivalently the dToda system 

Uti = -VT, Vt^ = -VUT, (2) 

and the 1-layer Benney (-B) system 

ut2=2{uut^+vt^), vt^ =2{uv)t^. (3) 

These systems represent themselves two distinguished examples of integrable 2-component hydro- 
dynamical type systems (see e.g. [5]). The dToda equation is the 1+1-dimensional version of the 
Boyer-Finley equation from the general relativity [6]. It also arises in various problems of fluid 
mechanics (see e.g. [7, 8, 9]). The B system describes long waves in shallow water with free sur- 
face in a gravitational field [10]. It represents the dispersionless limit of the celebrated nonlinear 
Schrodinger equation [11]. Recently, the B system became a crucial ingredient in the analysis of 
the universality of critical behaviour for nonlinear equations [12]. In general, the B system (3) is 
an excellent laboratory for study properties of integrable hydro dynamical type systems. 

Thus, the analysis of the singular sector for the H model in the large N limit reduces to the 
corresponding analysis for the dToda hierarchy. Recently the authors [13] have provided a method 
for studying the hodograph solutions and their singular sectors (gradient catastrophe points) of 
the family of coupled KdV hierarchies dcKdVm i'm > 1), where the case m = 2 corresponds to 
the B hierarchy. The study of these singular sectors was already initiated in [14]. In the method 
of [13] the hodograph solutions are described by the critical points of scalar functions that satisfy 
Euler-Poisson-Darboux (EPD) equations [15]. This property simplifies drastically the analysis and 
classification of the singular sectors . 

In the present paper we show that the results of [13] can be extended to the dToda hierarchy 
thus providing us with complete analysis of the singular sector of H model for the zero genus case. 
Moreover, we demonstrate that there is deep similarity between the properties of the sets of the 
hodograph solutions for the dToda and B hierarchies. 

Our main observation is that in both cases the hodograph solutions represent the critical points 

dW 



of a scalar function W which depends Hnearly on the coordinates t = (to, ti, ^2, • • •)) where t is 
the set of coupling constants for the H model (to = T) and the set of flow parameters in the case 
of the B hierarchy. Moreover, these functions satisfy the Euler-Poisson-Darboux equation E{a, a) 
[15] 

l iorWB 

^,JPw (dw dw\ I 2 

-- for Wh- 

V z 
Both functions are related according to a transformation [15] which maps solutions of £J(— 1/2, —1/2) 
into solutions -E(l/2, 1/2) 

where the parameters of both models are identified according to 

tB,n = {n + 2)tH,n+2, n > 0. 

The equation (5) and its multidimensional version are well known for a long time in the classical 
geometry [15]. Its relevance to the theory of Whitham equations has been demonstrated recently in 
the papers [16]- [17]. The observation made in [13] that the hodograph equations for the B system 
(3) have the form (4)-(5) seems to be new, though results close to this have been provided in the 
papers [12, 18]. 

The paper is organized as follows. In Section 2 we show how the hodograph equations of the 
B system and the H models are the equations for the critical points of certain solutions of Euler- 
Poisson-Darboux equations. In Section 3 we use the Euler-Poisson-Darboux equations to formulate 
a common description of the singular sectors for theses hodograph equations. Our results are 
summarized into three Propositions which exhibit the deep connection arising between the singular 
sectors of both models. Moreover, the corresponding singular classes are also characterized in terms 
of the behaviour near A = /3j (i = 1, 2) of the S*- function [B system) and the eigenvalue density p 
[H model). Section 4 is devoted to a method for the explicit determination of singular classes by 
means of constraints for the t-parameters. Finally, some illustrative examples are provided. 

2 Hodograph equations and the Euler-Poisson-Darboux 
equations 

B hierarchy 

The B system (3) is a member of a dispersionless integrable hierarchy which describe deformations 
of the curve (see e.g. [14, 19]). 

/ = (A-/3i)(A-/32). (6) 

The flows 

/3(t), t = {x := to,ti,t2, . . .) 

are characterized by the following condition: There exists a family of functions ^(A, t, /3) satisfying 

a5(A,t,/3(t)) 



dtn 
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fl„(A,/3(t)), n>0. (7) 



where 



^n(A,/3) 



i,n+l 



=) V(A-/3i)(A-/32). 



(8) 



V(A-/3i)(A-/32) 

Functions S" which satisfy (7) are referred to as action functions in the theory of dispersionless 
integrable systems (see e.g. [20]). Notice that for n = Eq.(7) reads 



so that the system (7) is equivalent to 



dS_ 
dx' 



^- d n 

Otn 



and, in terms of /3, it can be rewritten in the hydrodynamical form 



dtr, 



f^„(A,/3) 



A=ft 



d^Pi, i = l,2. 



(9) 



(10) 



The ti-flow of this hierarchy is the B system (3) (t := ti), which in terms of /3 reads 

9i/3i = ^(3/3i + /32)/?i., 

9t/32 = ^(3/32 + /3i)/32x. 

It was proved in [13] that the system (4) for the function 

dA AVB(A,t) 



(11) 



WB{t,P) := ^ 



where 



7^^^ V(A-/3i)(A-/32)' 
VB(A,t) = J]A"t„, 



(12) 



n>0 



is a system of hodograph equations for the B hierarchy. Moreover, the action function for the 
corresponding solutions is given by 



n>0 



where 



S{X, t,(3) = Y^ in 5^„(A, (3) = hniX, t, (3) V(A - /3i)(A - /32). 

XVB{X,t) 



(13) 



hB{X,t,f3):-- 



V(A-/3i)(A-/32) 



Obviously, Wb satisfies the Euler-Poisson-Darboux equation E{^, ^). Written explicitly, Wb rep- 
resents itself the series 



Wb = ^(A + /32) + I (3/?? + 2/3i/?2 + 3/3|) + ^ (5/3? + 3/3^/32 + 3/3i/3| + 5/3|) 



X 

2' 



+ ^(35/3f + 20/33/32 + 18/32/31 + 20/3i/3| + 35/3|) + ■ 



(14) 



The hodograph equations (4) with t„ = for n > 4 take the form 

8x + 4ti(3/3i + /32) + 3*2 {5(3l + 2/3i/32 + /3|) + | (140/3? + 60/3?/32 + 36/3i/3| + 20/3i) = 0, 

o 

8x + 4ti(/3i + 3/32) + 3t2 (/3f + 2/3i/32 + 5/3|) + | (140/31 + 60/?|/?i + 36/32/?? + 20/3?) = 0. 

(15) 
Detailed analysis of equations (15) will be performed in section 4. 



dToda hierarchy 

It is known (see e.g. [1]) that as A^ — )■ oo the asymptotic density of eigenvalues {Ai, . . . , Xn} of the 
H model in the zero genus case concentrates on a single interval [/3i,/32] and is given by [1] 

p(A, t, P) := ^^^ V(A-/?i)(A-/32), 
where 



V(A-/3i)(A-/32)' 

Here © denotes the projection on strictly positive powers of A. Moreover, the endpoints /3i and /?2 
of the eigenvalue support are determined by the equations 

dA dMX,t) _^^ /dA XdMX,t) _ ^^^^ 



f^2i7T ^(A-/3i)(A-/32) Jy^in ^(A - /3i) (A - /32) 

where 



Fj^(A,t):= J^A"t„, 



n>l 



and 7 denotes a large positively oriented circle |A| = R. It is immediate to see that these equations 
are equivalent to the system (4) for the critical points of the function 

WHit,(3) := T (/3i + (32)+<l^ ^{X - /3i) (A - ^2) dxVniX, t), (17) 

which satisfies the Euler-Poisson-Darboux equation E^—-^, — ^)- It is also easy to check that the 
functions Wb and Wh are related according to 

Wb=4: .^ " -tH,l, tB,n = {n + 2)tH,n+2, U > 0. (18) 

OPl Op2 

The first terms of the function Wh are 

Wh = T(/3i + /32) - |(/3i - 132? ^ |(/3i + /32)(/3i - 132? ~ ^ (5/3? + 6/32/3i + 5/3|) (/3i - 132? 
- ^ (7/3? + 9/32/3? + 9/3i/3i + 7/3|) (/3i - /32)2 + • • • , 



The corresponding equations (4) for t„ = 0, Vn > 4 take the form 

(5/33 - 3/32/3? - /3i/3i - Pi) = 0, 



T + |(/32 - /3i) + I (-3/32 + 2/32/3i + /3|) 



t2 



3*3 



J + |(;3i - /32) + I (/3? + 2/32/31 - 3/32) + ^ (/33 + (j^m + 3/3|/3i - bpl) = 0. 

It is known that the hodograph system (16) provides solutions of the dispersionless Toda hier- 
archy (see e.g. [3, 4]). Indeed, for "times" T and ti one gets 



5ti/3i 



1 



(/3i-/?2)/3: 



IT, 



(19) 



which in terms of variables 



9t,/32 = -(/3i-/32)/32r, 



^ = ^(/5i+/32), ^ = ^(/3i-/32)', 



(20) 



becomes the dToda system (2). Moreover, for times ti,t2 the hodograph equations (16) imply the 
Benney system (3). If we use the dependent variables {u,v) and introduce the function 



f(t,u,v) := WH(t,u,v) -2Tu + 2tiv + At2UV, 
then the hodograph equations read 

(2T-4vt2 + duf = 0, 
\-2ti-4ut2 + d^f = 0, 



(21) 



(22) 



which, under trivial rescalings and the substitution {u,v) 1— )■ {v,u), coincides with the hyperbolic 
version of hodograph system (2.11) of [12]. In particular, from the EPD equation (5) it follows 
easily that 



d', 



d^ 



gp2 



lU^ 



di 



0. 



3 Characterization of singular sectors 

Using (4) and (5) we may now analyze the structure of singular sectors of the B system and the 
H model in a unified way. Thus, let us denote by Ai the set of solutions (t,/3) (/3i ^ /32) of the 
hodograph equations (4). There is a partition of Ai into a regular and a singular sector 



7W=7Wregu_y^smg^ 



such that given (f , P) £ A4 



(t,/3) E A^-g If det (^-^M^) / 0, (t,/3) G A^-ng ,f ^et f C^^^ ) = 0. 



9/3,9/3,- 



5/3,9/3, 



The elements of A^^^S, correspond to the case when the system (4) is uniquely solvable and hence, 
it defines a unique solution /3(t). The singular sector A^^^^S contains the degenerate critical points 
of the function W on which the implicit solutions (3{t) of the hodograph equations exhibit gradient 
catastrophe behaviour. 

The Euler-Poisson-Darboux equation is of great help to analyze the structure of A^^^'^S. Indeed, 
if (t,/3) G A^, as a consequence of (5) it is clear that 

0. (24) 



Consequently 






Thus, we have 



Proposition 1. Given {t,f3) G Jv[ then 

1. {t,/3) G A^^e^ ^f^^^ ^^ly ^flPw ^ Qand^ / 0. 

2. (t,/3) G Jvi^^^^ if and only at least one of the derivatives ^r-p-, a^ ^ vanishes. 

Furthermore, using (5) it follows at once that at any point (t, /3) G A^ all mixed derivatives 
5o d-'a W can be expressed in terms of linear combination of diagonal derivatives d"^ W and d^W. 

Let us now define A4ni,n2 as the set of points (t,/3) G A4 such that 

x^ / 0; r- = 0, Vl< /c < nj + 1, (i = 1, 2). (26) 

It is clear that 

-MS n-^nl°4 = f3, for (ni,n2) / {n',,n'^) 

Then, it follows at once that 

Proposition 2. The singular sector satisfies 

M''^9= J A^f;^5. (27) 

«l+n.2>l 

According to (18) and using (5) we have 

2 /OWh dWH\ 

/3i - /32 V dl3i 5/?2 / 

It is then immediate to conclude that 

Proposition 3. The singular sectors of the H model and the Benney system, with 
tB,n = (n + 2) tH,n+2 for n>0, satisfy 

(28) 
l-^S+i,n.+i C M'Zin, for m > l,n2 > 1. 

7 



Singular sectors in the H model 

In the applications of the H model in quantum gravity [1] the singular sectors are described in 
terms of the behaviour of the eigenvalue density p{\) at the endpoints (/3i,/32) of its support. We 
will now show how this description derives in a natural way from the classification (27) based on 

the subsets M.tt„ „, of the singular sector . To this end let us consider the derivative do^ Wh 
with A; > 2 , it is obviously proportional to the integral 



dA (,_fe) ^^.^/n/;^ - f.) (> - fe) ^ i ii - ft) '■»(^) , ^ ( (,_ft),„(,) 



2i7r ' '^' {\-PiY A2i^ (A-/3i)fc {k - l)\ 



sing 



A=/3i 



and a similar result follows for the derivatives dJ" W with k > 2. As a consequence we have 



Proposition 4. A point {t,f3) G Mh belongs to the singularity class -A4//„ „ if and only if 



2nj + l 



A=/3l 



/9(A,t,/3)~(A-ft)^ asA^A, (i = l,2) (29) 

This property is crucial to establish the relationship between the regularized singularity sectors of 
the H model and the minimal conformal models {p,q) with q = 2 [1]. 

Singular sectors in the B system 

In analogy with the H model we may characterize the classes -M^n n '^^ ^^^ singular sector of 
the B system in terms of the behaviour of S{X) at A = /3j (i = 1, 2) . Indeed the derivative dg'^ Wb 
with A; > 1 is proportional to the integral 

dA xVB/^{\-Pi){\-h) f dA {X-P2)hBiX) dl( ^ 

,^^^-^^^ (A - /3i)^+i = l2l^ (A-/3i)^+i =^[{>^-p2)hB{X)) 

and a similar result follows for the derivatives dg^ Wb with k > 2. As a consequence we have 
Proposition 5. A point {t,(3) G Mb belongs to the singularity class A^^„^„2 ^^ '^^^ '^''^^y ^-^ 

2n,+3 

5(A,t,/3)-(A-A)^ as A ^ ft, (i = l,2) (30) 

This result shows the existing duality between the eigenvalue density p of the H model and the 
S function of the B system. 

4 Explicit determination of singular sectors in terms 
of constrains fk{t) = 

Sill 2! 

The singular classes A4ni,n2 can be determined by means of systems of ni + 712 constraints for the 

coordinates t. To see this property notice that the points (t, (3) of A4ni,n2 are characterized by the 

equations 

d'^W 

--^ = 0, Vl<A;<n, + l, i = 1,2, (31) 



and 



d 



\ni+2 



w 



/O, i = l,2. 



Furthermore, the jacobian matrix of the the system of two equations 



is not singular as 



A := 



5/3f^+^ 


J, — ±,z. 


9"i+2 w 


Qm+2 ^^ 


9/3i"i+2 


9/3i9/3f+^ 


a"i+2 w 


5"2+2 1^ 


dp'^^+^dh 


9/?2"'+2 



(32) 



(33) 



/O. 



(34) 



Indeed, we notice that as a consequence of (5) the derivatives outside the diagonal of A are linear 
combinations of the derivatives {da, W, 1 < A; < n-j + 1, i = 1, 2}, so that from (31)-(32) we have 

Qni+2 ^r Qn2+2 y^ 

A = : • ; / 0. 



d^l^ 



ni+2 



9/3"2+2 



Therefore, one can solve (33) and get a solution /3(t). Substituting this solution in the remaining 
equations (31) gives n\ + n2 constraints of the form 

fi,{t) = 0, A; = l,...,ni +71,2. 

It is not difficult to determine the solutions of (31)-(32) when t„ = 0, n > 4. Lists of such solutions 
are given next. 



Example: the H model 

Let us consider the H model and assume that t„ = if n > 4 . Then it follows that the condition 
det(9/3-/3^. WH{t,(3)) = which characterizes M^ reads 

32*2 _ 24^i3 + 72t2ts{^i + 132) + 27tl {^f + 6/?2/3i + Z?!) = 0. 

There are two generic classes in A^ji^ given by 7Wj:^-i^q and TWji^q^. Assuming that (ti, t2,i3,/3i,/32) 
must be real, both classes include two cases 



■^f/,1 ■ { 



1) X 



/3i 



2tl - 9^/&T^4^^ 

6t3 



2t2 + ^sertg/^ ^ 

^^^ ' ^2 



-2t2 + 3v/36rto/^ 



6t3 



2) X 



Atl + g^^er^ (i - i^) 4^' 



12t3 



-4t2+ (l + i\/3) v^^ertf'^ -4i2-3(l + i^/3) v'^ertg 

Pi = TTTI > P2 = 



2/3 



12t3 
where we assume ^3 < in the second case. 



12t3 



M 



smg 
i/,0,1 



1) X 



/3i 



2) X 



2ti - 9v/6T2 t^^ 

-2t2 + 3v^^6Tt3^'^ 
6t^ 



, h 



2t2 + v^^6Tt3/^ 

&f3 ' 



f3i 



4ti + 9^/6T^(l-i^/3)t^/' 

12t3 

-4*2 -3(1 + iVs) v'^sertg/^ 

12t^ ' 



/32 



-4t2 + (i + i\/3) v'^ert 

12*;^ 



2/3 



where we assume ^3 < in the second case. 

Example: the B system 

Let us now consider the system of hodograph equations for the B system with f„ = for all n > 4. 
Now the condition det(9^^^,WB(*i/3)) = reduces to 

32*2 ^ 96^2 (/3i + P2) + 702*2/32/32 + 72 (3t^ + ht^) /3i/32 + 12 (3*1 + 13*1*3) {^f + /3|) + 
486*2*3 {Ml + /3i/3i) + 90*2*3 (/3? + /3f ) + 180*^ (/32/3? + /3|/3i) + 45*1 (/?! + /^l) = 0. 
There are two classes in 7W^ ^ q 
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1. X 



p 



-45*3*^ + 180*1*1*2 + \/l5(8*i*3 - 


-3*2)^*2(3*2-8*1*3) 




360*1 


■) 




5*2*3 + Vl^Jtl (3*i - 8*1*3) 
20*2 


-3*2*3 + Vl5\/tl (3*2 - 


- 8*1*3) 


-45*3*1 + 180*1*1*2 - \/l5(8*i*3 - 


-3*2)^*2(3*2-8*1*3) 




360*3 






-5*2*3 + \/l5^ii (3*i - 8*1*3) 
20*2 


3*2*3 + \/l5y'*i (3*i - 
^' " " 12ti 


- 8*1*3) 



2. X 



/3i 
In the same way it follows that A^^ g 1 ^^^ ^^° classes given by 



1. X 



P 



2. X 



h 



-45*3*1 + 180*1*1*2 - \/l5(8*i*3 - 


-3*1)^*2(3*1-8*1*3) 




360*1 






3*2*3 + V^\jtl (3*2 - 8*1*3) 
12*2 


-5*2*3 + V^yjtl (3*1 - 
^' 20*2 


- 8*1*3) 


-45*3*1 + 180*1*1*2 + \/l5(8*i*3 - 


-3*2)^*2(3*1-8*1*3) 




360*1 






-3*2*3 + \/l5 Y^*i (3*1 - 8*1*3) 
12*2 


5*2*3 + \/l5^*i (3*1 - 
^^ - - 20*2 


- 8*1*3) 
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